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Abstract 

After a more general assumption on the influence of the bulk on the brane, we extend some 

n 



conclusions by Maarlons tl at. [jj and Santos el ai 2j on the asyioplol ic behavior of Bianchi I 
brane worlds. As a consequence of the nonlocal anisotropic stresses induced by the bulk, in most 
of our models, the brane does not isotropize and the nonlocal energy does not vanish in the limit 
in which the mean radius goes to infinity. We have also found the intriguing possibility that the 
inflation due to the cosmological constant might be prevented by the interaction with the bulk. We 
show that the problem for the mean radius can be completely solved in our models, which include 
as particular cases those in the references above. 
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I. INTRODUCTION 



The main stream approach to try and solve the problems arising from the break down 
of Einstein's theory at high energies is that of considering it a particular limit of a more 
general theory, as happened with Newton's theory in relation to that of Einstein's. One 
of this schemes is the brane gravity picture 0, 0], according to which matter fields are 
confined to a hypersurface (the brane) with three spatial dimensions embedded in a higher 
dimensional space (the bulk) on which gravity can act also. 

Considerable efforts in this area of research have been directed to testing that scheme by 
deducing cosmological implications from it. The usual framework for these studies is the 

n 

geometric approach of 5], which provides the gravitational field equations induced on the 
brane along with conservation equations for bulk degrees of freedom. Reviews on this topic 
can be found, for instance, in 0, Q, 0| . 

Studies concerning FRW universes [l| hinted that the dynamics of the early Universe in 
this new scenario is rather peculiar because of the modifications in the Friedmann equation. 
Then, when anisotropic models were considered ^(|, it was shown that discrepancies with 
respect to the standard scenario extended to the shear dynamics as well. 

In the set-up of general relativity, Wald nj showed that a positive cosmological constant 
leads to the isotropization of expanding homogeneous cosmological models . The anisotropy 
dissipation in brane world inflation in the absence of effective cosmological constant was 
analyzed by Maartens, Sahni and Saini [J. Taking into account the cosmological constant, 
Santos, Vernizzi and Ferreira derived a set of sufficient conditions that allow extending 
Wald's result to the brane world scenario Q (see also Q for an alternative view on the 
same problem). The last three works were done directly on the brane, where an additional 
hypothesis has to be done for the nonlocal anisotropic stresses induced by the bulk, for they 
are not given by any evolution equation on the brane. This naturally raises the question of 
to what extent the results depend on the aforementioned hypothesis. 

The goal of this work is to extend the results in by considering a more general 

additional hypothesis. We will analyse, as done in those works, a Bianchi I brane embedded 
in a 5-dimensional bulk. We assume that the nonlocal stresses satisfy a condition which 
reduces to the one considered in Refs. in particular cases. We will show that the 

evolution equations on the brane can be integrated and conclude that, generically, we have 
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exponential inflation or an asymptotic power law for the mean radius, but the models do 
not isotropize unless they belong to the class considered in the aforementioned works. 



II. BIANCHI I BRANE MODELS 

Let us consider a Bianchi I brane with the induced metric 

ds 2 = -dt 2 + a\(t) dx 2 + aj(t) dy 2 + a\(t) dz 2 (1) 

when the matter on the brane is a perfect fluid of density p and pressure p. We define as 
usual the mean radius a = {a\a2a^) 1 ^ and the mean Hubble parameter H = a/a = ~ X)i Hi, 
with Hi = a,i/ai. 

We will be using the equations induced on the brane derived by Shiromizu, Maeda and 
Sasaki ISfl but follow the notation of . For the metric (|T|), the dynamics of the mean 
radius a in this model can be described by the following equations on the brane: 

3H + 3H 1 + „'< + -jLw = A - y (p + 3p) - ^(2p + 3p)^, (3) 

6 

&^ + 3Ha^ = —^rVfu,, (4) 
U + AHU + a^Vfu, = 0,, (5) 

where a dot denotes i^V^, with = d/dt. Throughout the paper k 2 /8tt and k 2 /8ir will, 
respectively, be the effective Newton constant on the brane and in the bulk ; A will be the 
effective 4D cosmological constant, and the tension on the brane will be A = 6k 2 /k 4 . On 
the other hand, U stands for the effective nonlocal energy density, and a^ v denotes the shear 
scalar, which satisfies (HJ) and 

( r 2 = a^a^ = j2(H l -H) 2 . (6) 

i=l 

In addition, if w M is the four velocity of an observer on the brane comoving with matter, 
then h^y = g^ u + u^u u projects into the comoving rest-space, and U and V^ v will be related 
to S^v, which is the projection on the brane of the 5D Weyl tensor, through 

U = -^^uV, (7) 
o 

TV = (hfhf - Sap. (8) 



There is also, in principle, the constraint that the projected spatial covariant derivative 
of Vfu/ vanishes: 

D*7V = 0, (9) 

but this result holds identically for the metric ©. On the other hand, as a consequence of 
system (0)-©, we get the conservation law 

p + 3H(p + p) = 0. (10) 

III. EXACT EXAMPLES AND ASYMPTOTIC BEHAVIOR 

Since there is no evolution equation on the brane for the nonlocal anisotropic stress V^, 
but only the constraint (j9*)l. which in this case does not provide any information, some 
additional hypothesis is necessary to integrate the system (J2J)-©- 

In order to get some insight on the problem, we will restrict ourselves for a moment to 
the vacuum case p = p = 0, where we get from (J2J)-(I3J) 

U ='f±(H + ? > H 2 -k). (11) 

Let us consider under which conditions the stable asymptotic behavior may approach a 
power law a ~ t k , {k > 0). In such a case H , H — > 0, U — > — /t 2 AA/2 and, because of (JHJ, 

o^Vp, -> 2k 2 XAH. (12) 

We will show in the following that, if one chooses the unknown V^ u so that this condition 
is satisfied, the asymptotic behavior described by a power law is stable. Furthermore, this 
will happen even in the presence of a fluid. 

So as to integrate the system (J2JH© one often assumes U = or the more general 
g^V^v = 0, which has been used in 2j to discuss the stability of the de Sitter spacetime. 
It has been proved Q| that spatial homogeneity follows from the integrability conditions 
for vanishing non-local anisotropic stress and energy flux, which are two of the three bulk 
degrees of freedom. Therefore, the considerable simplification arising from switching off 
those quantities is consistent with having a Bianchi I metric on the brane. 

It must be pointed out that other choices can be found in the literature; Barrow and 
Maartens in an investigation of early times shear anisotropy in an inhomogeneous 
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universe, assumed that behaves qualitatively like a general 4D anisotropic stress, in 
particular they chose to be proportional to the energy density of the anisotropic source 



(which we will denote with p) so that V^ v = pC^ v with C^ v = 0, \jC^ v C^ v = 0(1), and 

Although the asymptotic behavior (j!2)) may happen with many choices of cr^V^, we 
will consider only a simple family of models. In the following, we will assume that u^V^ 
is proportional to the Hubble parameter, so that it can be written as 

o^V^ = 2k 2 \T(cl) H, (13) 



m 



[JUL11U 



which reduces to the cases discussed in |l|,|2| when the function T(a) vanishes. By making this 
hypothesis, we use a less restrictive assumption to check whether the conclusions reached 
hold in a more general context, while still being able to integrate the evolution 
equations. Note that |T2|) is recovered, even before the asymptotic behavior is reached, 
provided T(a) = A. 

Exactly as happens with the restrictive T = case often used in the bibliography, it 
remains an open problem whether the more general assumption we make is compatible with 
a full 5D solution (see for a recent attempt at tackling the 5D problem). Judging on the 
suitability of such guesses is, therefore, not possible for the time being, but they definitely 
help unveiling possible unexpected consequences of the interference between bulk and brane. 

The nonlocal energy density U is given by the conservation equation (JHJ), which is easily 
solved under the assumption (JT3J): 

U = — ^ - J a 3 T{a) da, (14) 

where uq is an arbitrary integration constant. We see here that, for expanding universes, 

lim U = ——n 2 X lim T(a), (15) 

a— *oo 2 a^roo 

(provided the limit on the right side exists or is infinite) so that it does not always vanish 
asymptotically, although it will go to zero in the special cases in which T(a) — > (which 
of course include the choice T(a) = of Refs. [l|, Q), where it reduces to the result by 
Toporensky Expression (fTB} suggests there is the possibility that, as a consequence 

of the interaction of the brane with the bulk through V^, the asymptotic vanishing of U 
arising in the cases studied in Q might be evitable. 



By contracting with a^ v and using (jl3j) . one readily gets 



for any constant <To, so that 



a 2 = ^f + ^ I a 5 T(a)da, (16) 



lim a 2 = 4 lim T(a). (17) 

a— >oo a^oo 



Then, it is clear that asymptotically V must be non-negative. If T — * the brane world 
isotropizes as a — > oo, as described in , but, for all other asymptotic behaviours of T, 
there will be a remaining anisotropy induced by the nonlocal anisotropic stresses. Note that 
there is also the possibility that the anisotropy grows with a, which would correspond to 
U < at late times. Models which do not isotropize in a recollapsing situation, unlike ours, 



were found in [17(, and they were characterized by V^ u and U < 0. 

Clearly, our work and many others suggest that the interplay between the bulk degrees 
of freedom and the dynamics on the brane is rather non-trivial. An interesting illustration 
of this is a recent work |Tj| which investigated the dynamics of a flat isotropic brane world 
with a perfect fluid with equation of state p = (7 — l)p and a scalar field with a power-law 
potential. There it was found that the number and the stability of fixed points of the system 
describing the dynamics of the model would depend not only on whether U vanishes or not, 
but also on its sign. 

If we assume the equation of state p = (7 — l)p, the conservation law (fTT)|) is equivalent 

to 

' = (18) 

with an arbitrary constant po- 

If we insert (JHJ), (fTB^) and (JTSJ) in (J2J), we get the generalized Friedmann equation 

3H 2 = ^ + Po + _^L_ + f\ a fa 3 T(a)da] da. (19) 

a 4 2a 6 a 3 T 2k Aa 7 a 6 J [ J K ' \ K ' 



The orbits (H(U), cr(U)) were found in |2j for the special case T = 0, but we can see that, 
for any T(a), the problem for a(t), a(t) and U(t) may be completely solved from equation 
()19|) by means of a quadrature and a function inversion. The integrals may be explicitly 
computed (in terms of elementary or elliptic functions) for different choices for the constants 
Mo, Cq, po, 7 and A and the function T(a). The simplest cases are solutions with no fluid 
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(po = 0). If r(a) = A, by using the parameter < w < oo one may write the solution as 



48a/2m 3 



jL= (sinh w-w), (20) 



o 



a ° shm™, (21) 



2v^ 2 ' 

which at late time corresponds to a ~ t 1//2 . 

In contrast, if «o = po — and T(a) = a for a constant a < A, we get 

a 6 = ^ 

2(A - a) 

Other exact solutions with uq = can obtained with dust (7 = 1) or a stiff fluid (7 = 2). 
Let us look first at the dust cases, which, arguably, are very interesting from the observational 
point of view Q|. For T(a) = a we have 



sinh 2 J3(A - a) t. (22) 



, Po . u2 y3(A-a)t 

a = sinh h 



A - a 2 \ 2k 2 X{A - a) 



Po + K 2 \al 



sinhA/3(A-a)t. (23) 



20. 211 that has not been 



The latter is a generalization of the Heckmann- Shucking metric 
discussed in the literature so far. Clearly, the a — > A limit of the solution ()23|) is regular, 
with the form of a second-order polynomial in t. The effects of the shear and the quadratic 
corrections are of the same order, they dominate at early times, and a ~ t 1 ^ 3 . In contrast, 
at late times the model is neither aware of the anisotropy nor of the extra dimensional 
ingredients, and log a oc t, or a ~ t 2 ^ 3 if A = a. On the other hand, for a stiff fluid (7 = 2), 
and r(a) = a, we have 

„6_ 2po + ^ 2 • ,2 KT K —+ , Po 



. sinh 2 J3(A - a) t + -= r sinh 2J3(A - a) t. (24) 

2(A-a) V 1 J ^/2^A(A-a) V 1 ^ 1 ; 

This solution too has a regular a — > A limit in the form of a second-order polynomial in t. 
Like in the dust solution above, the effects of the shear and the quadratic corrections are 
equally important, but they dominate at late times, instead; in that regime we either have 
log a oc t, or a ~ t 1 / 3 if A = a. On the contrary, at early times the model becomes isotropic 
and relativistic, and a ~ t 1 / 6 . 

We can also obtain directly from (JT§j) some general late-time results of interest. If asymp- 
totically r(a) ~ a is constant, then for A — a > 0, the brane inflates exponentially when 
a — > 00, as happened for a = 0, A > 2}. But now there is another interesting possibility: 
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the repulsion due to the cosmological constant may be neutralized by the bulk influence 
in which case the asymptotic behavior of a is a power law. For simplicity we show this in 
vacuum and write in terms of the variable \x = —2H/3H 2 , which becomes a constant 
Hq for the power-law solutions a oc t 2 / 3 ^ , so 



fi + 



3Hfi- 



8HU + 2a^P ia , 



(2-//)=0. (25) 



2U + k 2 \A 

The existence of the constant solution, /i = /x , in the last equation requires that 



a^P^ = 3H 



"( U+ — )-3 U 



(26) 



this means that, for any solution of (11) which asymptotically approaches to a power law, 
then the quantity a^ u P^ v asymptotically behaves as Now, inserting (|2Tj) in (J3J) and 
(HJl, we get the asymptotic expression for the effective nonlocal energy density 

n 2 \a 1 9 , , . 
W = —= k AA, 27 

and the shear scalar 

, <x 2 12^(^o-4/3) 
^ (2-p )a 3 - + 4A > ( 28 ) 

where a and a\ are arbitrary constants. In addition, introducing (|2l)j) into it reads as 



fi = -3H(iM-iM )(2-fx). (29) 

It is easy to see that, for ordinary fluids, /i < 2, the solutions of Eq. (j2*U|) converge, respec- 
tively, to the stable fixed point /i = /io, which describes the asymptotic power-law solutions 
a = t 2 / 3 ^ . For instance, if we additionally assume a perfect fluid source, after using (|18jh 
()27|) and (|2"H|). the Einstein equation (J2J) becomes 

2 = 6a erg 4a p p 2 

a 3 ^o ^ 2a 6 (2 - fi )a^ 2k 2 \a^' 1 ' 

Hence, for 7 > po, asymptotically the brane will remain anisotropic (a 2 ~ 4A) and expand 

as a ~ t 2//3At ° or expand as a ~ t 2 / 37 for 7 < po- Note that, the exact solution given by 

(J2UJ) and (|2~Tj) for a model with no fluid, is also the approximated late time solution when 
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/ uo = 4/3<7 or 7 = 4/3< / uo- in both cases asymptotically the scale factor behaves a 
a ~ t 1 / 2 . 

In we assumed that a^ u P^ v depends on a and U. We could instead assume that it 
depends only on a so that asymptotically (J 13)) holds with 

r(a)=A + A, (31) 

for some constants /3 and 5 > 0. It is easy to see that in this case also the asymptotic 
behavior of a is a power law. For instance, if we additionally assume that 7 > 4/3 (or that 
there is no fluid), asymptotically the brane will remain anisotropic (a 2 ~ 4A) and expand 
as a ~ 

t 2/38 for 5 < 4/3 an d as a ~ t 1 / 2 for 5 > 4/3 (or (3 = 0). For 7 < 4/3, one would have 
a ~ t 2 / 3u , with z/ = min(7, 5). 



IV. CONCLUSIONS 

By making the more general assumption (|T3*|) on the (unknown) influence of the bulk on 
the brane, we have shown that some conclusions on the asymptotic behavior of Bianchi I 
brane worlds in Refs. can be generalized. Due to the nonlocal stresses, in most of our 

models, the nonlocal energy does not vanish in the limit a — > 00, and the brane does not 
isotropize. We have also found that, although nearly all our models inflate, there also exist 
the possibility that the inflation due to the cosmological constant might be prevented by the 
interaction with the bulk. Finally, we have shown that the problem for the mean radius a 
(as well as for a 2 and IX) can be completely solved in our models, which include as particular 
cases the Bianchi I branes for which the orbits and stability were analyzed in 
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